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Lee, Abramowicz & Kluzniak H2004fl demonstrated nu- 
merically that rotation of neutron star couples with os- 
cillations of its accretion disk, and excites resonances. No 
specific coupling was assumed, but magnetic field was sug- 
gested as the most likely one. Following this idea, we show 
(Petri YODSa paper I) that if the neutron star is non- 
axially symmetric and rotating, its gravity may provide 
the coupling and excite resonances. Here, we return to 
the original suggestion that the coupling is of a magnetic 
origin, and demonstrate how does it work in terms of a 
simple, analytic toy-model. 
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1. INTRODUCTION 

In some neutron star sources, the observed QPO frequen- 
cies obviously depend on the neutron star spin. For exam- 
ple, difference in frequencies of the double peaked QPO in 
the millisecond pulsar SAX J1808. 4-3658 is clearly equal 
to half of the pulsar spin (Wijnands et al. I2003|l . This 
made Kluzniak et al. I|2UU4|I to suggest that a resonance is 
excited by coupling accretion disk oscillation modes to the 
neutron star spin. The suggestion was fully confirmed by 
numerical simulations of the coupling (Lee, Abramowicz 
& Kluzniak, '2004') . It was found that a resonant response 
occurs when the difference between frequencies of the two 
modes equals to one-half of the spin frequency (as ob- 
served in SAX J1808. 4-3658 and other "fast rotators"), 
and when it equals to the spin frequency (as observed in 
"slow rotators" like XTE J1807-294). 

Lee et al. (|2004|l suggested that the coupling is most 
likely a magnetic one, but did not specified any concrete 
mechanism in their numerical simulations, introducing 
the coupling by a purely formal ansatz. Following their 
idea, we discussed in terms of a simple analytic toy-model 
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fPctri I^OSal Paper I) that also (non-axially symmetric) 
gravitational field of rotating neutron star may excite res- 
onances in accretion disk oscillations. 

Another point of view to account for the dichotomy 
between slow and fast rotators is given by Lamb & 
Miller l)2003|l who explain it in the framework of the sonic- 
point beat frequency model. 

In this Research Note, we discuss another toy model 
that provides the coupling by the neutron star magnetic 
field. We use the same notation as in Paper I. 

2. THE MODEL 

In this section, we describe the main features of the 
model, starting with a simple treatment of the accretion 
disk, assumed to be made of non interacting charged sin- 
gle particles orbiting in the equatorial plane of the star. 
Magnetohydrodynamical aspect of the disk such as pres- 
sure and current are therefore neglected. Particles evolve 
in a perfectly spherically symmetric gravitational poten- 
tial. The asymmetry arises from an eccentric misaligned 
dipolar magnetic field corotating with the neutron star. 

2.1. Eccentric stellar magnetic field 

The periodically varying epicyclic frequencies are intro- 
duced by adding a rotating asymmetric dipolar magnetic 
field to the background gravity. Generally, when dealing 
with an oblique rotator, the location of the magnetic mo- 
ment fi generating the dipolar magnetic field coincides 
with the centre of the neutron star (supposed to be a per- 
fect sphere). In this paper, we lift this assumption and 
shift the location of the magnetic moment to a point 
rs{t) = {vs ^ 0,ips ~ fl^,t,Zs) inside the star such that 
||rs(t)|| = \/rf+^ < i?, where is the stellar ra- 
dius. We use cylindrical coordinates denoted by {r,(p,z). 
Nevertheless, the origin of the coordinate system coincides 
with the centre of the neutron star. Furthermore, in order 
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to compute analytically such kind of magnetic field struc- 
ture, we assume that the star is made of an homogeneous 
and isotropic matter everywhere with total mass and 
spinning around its centre at an angular rate = e^, 
aligned with the z-axis. The magnetic field induced by the 
dipolar source is therefore : 



B{r,Lp,z,t) = 

4 TT 



3{fi- R)R n 

^5 



(1) 



The vector joining the source point rs{t) = (rs,ips,Zs) to 
the observer point r — (?', (p, z) is : 



R{t) 



rs{t) 



(2) 



Using the cylindrical frame of reference, the distance be- 
tween source point and observer is : 



r'^ + Tg — 2 r Ts cos tjj + {z ~ Zg)^ 



(3) 



where the azimuth in the corotating frame \s ip = — 
17* t. The magnetic moment anchored in the neutron star, 
rotates at the stellar speed such that : 

/x(t) = ^ [sin X {cos (17* t) + sin (il* i) ey } -f cos x e^] (4) 

where the obliquity, i.e. the angle between /x and SI*, is 
denoted by x- Moreover, each component of the magnetic 
field can be expressed as follows : 
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fi cos X 



fj, ■ R — ^ [sin X (»" cos ip — rg) + cos x ~ ^s)] 



(6) 

(7) 
(8) 



The total linear response of the disk is then the sum 
of each perturbation corresponding to one particular az- 
imuthal mode m. Because the perturber is inside the star 
and the disk never reaches the stellar surface, the Fourier 
coefficients of each component B™ {i = i,(p,z) never di- 
verge. It is convenient to introduce the Fourier decompo- 
sition of the magnetic field component Eq. ©-Q order 
to describe the response of the test particle. Moreover, be- 
cause the evaluation of the Fourier coefficients requires to 
integrate terms containing cos (m ijj) in the integrand, the 
value of these coefficients decreases rapidly with increas- 
ing azimuthal number m. As a result, only low azimuthal 
modes will influence significantly the evolution of the disk. 
Keeping only the few first terms in the expansion is suf- 
ficient to achieve reasonable accuracy. For discussing the 
results, we only keep the three first modes, namely, the 
dipolar, quadrupolar and octupolar moments (m = 1,2,3 
respectively) . 



2.2. Equation of motion for a charged test particle 

All particles evolve in the gravitatomagnetic field imposed 
by the rotating neutron star. To keep things as simple as 
possible, their motion is described in the guiding centre 
approximation. The drift arising from the gyration around 
the local magnetic field is not an essential feature we want 
to discuss here. As a consequence, magnetic curvature and 
gradient as well as gravity drift motions are ignored in this 
study. Nevertheless, the main characteristic consisting of a 
periodic variation in the epicyclic frequencies is preserved. 
The equation of motion then reads : 



G = g 



GAB 



(9) 



where G is the location of the guiding centre and the dot 
means time derivative d/dt. The gravitational field of the 
star Af, is denoted hy g = V(GM,/\A^^"T^). The mass 
and the charge of the particle are denoted respectively by 
rric and Qc. If for instance the magnetic gradient drift is 
taken into account, the term — {fi^ / 2 nic B) 'V B^ should 
be added to the right hand side of Eq. © where fic is an 
adiabatic invariant, namely the magnetic moment of the 
test particle gyrating along the local field line. This would 
introduce another modulation of the epicyclic frequencies 
which is already included in the Lorentz force q^G f\ B. 
Thus, the physical behaviour is not changed by neglecting 
the drift motion of the guiding centre G. Expressed in 
cylindrical coordinates, Eq. I© develops into 



r — r ip 



.gr H {r ip Bz - z B^ 



2r ip + r (f = — ^ (i i?r ~ t' B^) 



5z H [rB^ -npBy) 



(10) 

(11) 
(12) 



The magnetic field B is assumed to be weak enough for the 
fiow to remain essentially hydrodynamical (weakly mag- 
netized thin disk approximation). We therefore treat B 
as a perturbation of order e <C 1. The perturbations in- 
duced in the fiow are of the same order of magnitude than 
B, i.e. of order e. The perturbed orbit and velocity of 
the test particle in the radial and vertical direction are 
also of order e, {r,z} = 0(e) whereas the azimuth varies 
as 93 « i7k ~ \/ G Af * / Tq + 0(e), the Keplerian orbital 
frequency at the radius of the orbit vq. In the equations 
of motion (fTn|l . (fTT|l and ifT^ . terms such as {r, i} x Bi 
with i = {i,p,z} are second order 0(6^) and we neglect 
them. According to this simplification, the right hand side 
of Eq. Hll|) vanishes. Eq. (|ll|l states the conservation of 
angular momentum of the particle and integrates into 
L = mr'^ ifi = const where L is the angular momentum 
of the particle. This is an obvious integral of motion for 
this problem (to the aforementioned approximation). 

However, we are only interested in the vertical motion 
experienced by the test particles in response to the per- 
turbed magnetic field. Indeed, the response of an accretion 
disc to an inclined rotating magnetic dipole has been stud- 
ied by Terquem & Papaloizou They showed that 
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the vertical displacement (i.e. the warping) is the domi- 
nant effect in the thin disc while the horizontal pertur- 
bations can be neglected. Perturbing H12|) and developing 
to first order in the perturbation around the equilibrium 
Keplerian orbit (contained in the equatorial plane) defined 
by (ro, ipo = fik t, zq = 0), the vertical motion reads : 

9o 



5z 



■ To f2k Br 



(13) 



To avoid mathematical irrelevant complications, we as- 
sume the rotator to be "aligned" with the star in the sense 
that fi and fl^. are parallel (x = 0). Therefore the radial 
component of the magnetic field reads : 



3 ^0 - Zs) {r - rs cos ip) 



At: 

By developing in a Fourier series, we obtain 



j-oo 



Br = (z - Zs) > B,"(r, z) cos (m ?/>) 



=0 



(14) 



(15) 



where the Fourier coefficients are given by 



i?r(r. 
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COS {mip) (iV'(16) 



periodically varying eigenfrequency which is also excited 
by a driven force. It is well known that some resonances 
will therefore occurs in this system. Namely, we expect 
three kind of resonances corresponding to : 

— a corotation resonance at the radius where the angular 
velocity of the test particle equals the rotation speed 
of the magnetic structure (which is the equal to the 
stellar rotation rate). Corotation is only possible for 
prograde motion. The resonance condition determining 
the corotating radius is simply ilk = 0*; 

— a driven resonance at the radius where the verti- 
cal epicyclic frequency equals the frequency of each 
mode of the magnetic perturbation as seen in the 
locally corotating frame. The resonance condition is 

m \^^, — fikl = Kzl 

— a parametric resonance related to the time- varying ver- 
tical epicyclic frequency, (Hill equation). The rotation 
of the magnetosphere induces a sinusoidally variation 
of the vertical epicyclic frequency leading to the well 
known Mathieu's equation for a given azimuthal mode 
m. The resonance condition is derived as followed : 



where is the Kronecker symbol. Note that B™ do not 
have the dimension of a magnetic field because of their 
definition Eq. (|15|l . These coefficients, which are function 
of the space position (r, z) decrease with increasing az- 
imuthal number m (for (r, z) fixed). Keeping the first few 
coefficients is sufficient to achieve a reasonable accuracy 
(we retain the first three). Putting the expansion Eq. (|15|l 
into the vertical equation of motion Eq. (|13|l . we get the 
fundamental equation to describe vertical forced oscilla- 
tions of a test particle as follows : 



(18) 



ni 



gc 

Too 



+ 00 



r!k J2 ^r(^o,^o) cos{m(r!k-f^*)0 



-z = roiik^s 

TOe 



-t-CXD 

E 

m=0 



B™(ro,2o) cos{m(r2k 



(17) 



Note that the Fourier coefficients i?™(ro,zo) in this last 
equation are evaluated at the location of the unperturbed 
orbit and do no longer depend on the perturbed position 
(r, z). To the lowest order of the expansion, this approxi- 
mation is justified. We recognize a Hill equation (periodic 
variation of the eigenfrequency of the system on the left 
hand side) with a periodic driving force (on the right hand 
side). 

2.3. Resonance conditions 

This equation is very similar to the one obtained in the 
case where solely gravity perturbation exists. The discus- 
sion is therefore exactly the same as in paper I. Here we 
recall the main results, adapted to the magnetic config- 
uration. Eq. l(T7| describes an harmonic oscillator with 



where n > 1 is a natural integer. Note that the driven 
resonance is a special case of the the parametric reso- 
nance for n — 2. However, their growth rate differ by 
the timescale of the amplitude magnification. Driving 
causes a linear growth in time while parametric resonance 
causes an exponential growth. We also rewrite the vertical 
epicyclic frequency as instead of fJk in order to apply 
the results to a more general case which could include 
stronger magnetic fields or general relativistic effects. 

Consequently, the resonance conditions for the mag- 
netized rotator are exactly the same as for the unmagne- 
tized rotator of paper I, as long as the oscillations remain 
in the linear regime (i.e. in the thin disk approximation 
for which the vertical motion of the particle remains small 
with respect to the radius of the unperturbed orbit). For 
more details and a discussion on these results, we refer the 
reader to paper I. 

The dichotomy between the QPOs in fast and slow ro- 
tators has been explained by Lee, Abramowicz & Kluzniak 
(2004) as a consequence of the fact that the coupling be- 
tween the neutron star spin and the modes of accretion 
disk oscillations excites possible resonances at different 
locations in the disk, either very close to the star, or far 
away from it. While we agree with this results, we point 
out that (as found in Paper I) the spectrum of modes that 
could be in a resonance is more complex than that con- 
sider by Lee et al. Their discussion was concentrated on 
the 3:2 resonance that occurs between the radial and ver- 
tical epicyclic modes only in strong gravity. We identified 
possible forced resonances between different modes, that 
may occur in both strong and weak gravity. Our analysis 
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was done in a linear regime, so, at the moment, we are only 
able to say that such resonances may in principle exist^. 

3. CONCLUSION 

The toy-model discussed in this Research Note illustrates 
the idea (Kluzniak et al., IMHI Lee et al. .2004 Petri 
l2UU5af that the double peak QPOs in neutron stars 
sources may be due to coupling rotation of neutron star 
sources to modes in accretion disk oscillations, and excit- 
ing resonances. 

The model is physically very specific. It explicitly 
shows how the rotating magnetic field of the neutron 
stars could couple with the disk dynamics and oscilla- 
tions. It confirms general discussion and numerical results 
of Kluzniak et al., I)2()()4|l and Lee et al. I|2()()4|l . in par- 
ticular the important one that the strongest resonant re- 
sponse occurs when the difference between frequencies of 
the two modes equals to one-half of the spin frequency (as 
observed in SAX J1808.4-3658 and other "fast rotators"), 
and when it equals to the spin frequency (as observed in 
"slow rotators" hke XTE J1807-294). 

When the MHD nature of the flow is taken into ac- 
count, QPOs can be explained by a mechanism similar 
to those exposed here fPetri [2005b|l . However, in an ac- 
creting system in which the neutron star is an oblique 
rotator, we expect a perturbation in the magnetic field to 
the same order of magnitude than the unperturbed one. 
Therefore, the linear analysis developed in this paper has 
to be extended to oscillations having non negligible am- 
plitude compared to the stationary state. Nonlinear oscil- 
lations therefore arise naturally in the magnetized accre- 
tion disk. Abramowicz et al. I|2003il showed that the non- 
linear resonance for the geodesic motion of a test particle 
can lead to the 3:2 ratio for the two main resonances. 
Nevertheless, an extension to 3D MHD fiow in curved 
spacetime is required to make quantitative accurate pre- 
dictions of the peak frequencies variation correlated with 
the accretion rate for instance. 
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